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Abstract. We present results concerning the linear (radial and non-radial) oscillations of non-rotating superfluid 
neutron stars in Newtonian physics. We use a simple two-fluid model to describe the superfluid neutron star, 
where one fluid consists of the superfluid neutrons, while the second fluid contains all the comoving constituents 
(protons, electrons). The two fluids are assumed to be "free" in the sense of absence of vortex-mediated forces 
like mutual friction or pinning, but they can be coupled by the equation of state, in particular by entrainment. We 
calculate numerically the eigen-frequencies and -modes of adiabatic oscillations, neglecting beta-reactions that 
would lead to dissipation. We find a doubling of all acoustic-type modes (f-modes, p-modes), and confirm the 
absence of g-modes in these superfiuid models. We show analytically and numerically that only in the case of 
non-stratified background models (i.e. with no composition gradient) can these doublets of acoustic modes be 
separated into two distinct families, which are characterised by either co- or counter-moving fluids respectively, 
and which are sometimes referred to as "ordinary"- and "superfluid" modes. In the general, stratifled case, 
however, this separation is not possible, and these acoustic modes can not be classified as being either purely 
"ordinary" or "superfluid" . We show how the properties of the two-fluid modes change as functions of the coupling 
by entrainment. We flnd avoided mode-crossings for the stratifled models, while the crossings are not avoided in 
the non-stratified, separable case. The oscillations of normal-fluid neutron stars are recovered as a special case 
simply by locking the two fluids together. In this effective one-fluid case we flnd the usual singlet f- and p-modes, 
and we also flnd the expected g-modes of stratifled neutron star models. 
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1. Introduction 

The study of stellar oscillations has proved very fruitful in 
improving our understanding of the inner structure and 
dynamics of stars (the terms helio- and astro-seismology 
have been coined), for which the oscillation modes can 
often be observed rather directly. The best developed ex- 
ample of this probing of the internal structure of an astro- 
physical body via its oscillations is probably the Earth. In 
the case of neutron stars, the observation of oscillations is 
unfortunately not possible in such a direct way, and has 
not yet been achieved. In practically all cases we can only 
observe the regular radio-pulses of neutron stars, which 
are virtually unaffected by its oscillations and give infor- 
mation mostly about their rotation rate. Nevertheless this 
field bears great potential interest: on one hand the better 
understanding of neutron star oscillations could eventu- 
ally help to elucidate the phenomenon of glitches, which is 
probably the most striking and puzzling aspect of observed 
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neutron star dynamics. This phenomenon still represents 
somewhat of a mystery, even though th e crucial role of su - 
perfluidity seems w ell established (see Link et al. (2000 ); 
[Carter et al. (2000 ) for recent discussions). On the other 
hand, several highly sensitive gravitational wave detectors 
are expected to reach their full sensitivity within the next 
few years, and neutron star oscillations represent one of 
the potentially most interesting sources of gravitational 
waves. Gravitational wave detection could open a new and 
complementary observational window onto neutron stars, 
which would allow us to learn much more about their in- 
ner structure and dynamics than it is currently possible 
with the purely electro-magnetic observations. 



Currently most studies of neutron star oscillations are 
based on simple perfect fluid models, which neglects the 
crucial importance of superfluidity in neutron stars. The 
presence of substantial amounts of superfluid matter in 
neutron stars is backed by a number of theoretical calcu- 
lations of the state of matter at these extreme densities 



e.g. see [Baldo et al. (1992[ ); pjoberg (1976D ), and by the 
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qualitative success of superfluid models to accommodate 
observed features of glitches and their relaxation. 

The first study to point out the importance of su- 
perfluidity for the oscillation properties of neutron stars 
was by Epstein (1988), who has argued in a local (sound 
wave) analysis that superfluidity should lead to new modes 
and modify the previously known modes. |Lindblom & 



Mende p" (1994| ) have argued further for the existence of 



these modes, but failed to find them numerically. Lee 
(1995|)]presented the first numerical results indicating the 
presence of new modes that did not exist in perfect fluid 
models, and the absence of g-modes which would have 
be en present in the non-super fiuid case. A local analysis 
by Andersson fc Comer (2001a ) has given further analytic 
evidence for the absence of g-modes in simple superfluid 



models. The relativistic numerical analysis by Comer et al 
(1999|)]has shown an effective doubling of acoustic modes 
in superfluid models with respect to the normal fluid case. 



Recently this work has been extended by Andersson et al. 
(2002|]]to include entrainment, and they have shown that 
avoided mode crossings occur when one varies the en- 
trainment parameter. The relevance of superfluid oscilla- 
tions for gravitational wave detection has been discussed 
by Andersson & Comer (2001b). Some studies have also 
started to look at oscillations of rotating superfluid neu- 
tron stars ([Lindblom fc Mendell 2000| ) and ^edrakian fc| 
Wassci|man (2000D. 



Despite the number of studies on oscillations of non- 
rotating superfluid neutrons stars, we think that this prob- 
lem still deserves attention and that several points needed 
to be clarified. In particular it is worth emphasising the 
importance of stratification for the nature of superfiuid 
oscillations, a point that has not yet been fully appre- 
ciated. We demonstrate here that only in non-stratified 
models can the eigenmode spectrum be separated into 
two families of modes, one of which is identical to the 
case of a normal-fluid star, while the other is charac- 
terised by counter-motion of the two fluids and vanish- 
ing gravitational perturbation. These two distinct fami- 
lies are usually referred to as "ordinary" and "superfluid" 
modes. Stratification of the background star, however, 
couples these distinct mode-families and renders them 
non-separable. As a consequence every mode shares qual- 
itative properties of both families to some extent, and the 
resulting mode spectrum consists of modes that bear no 
direct connection to the normal fiuid case. The two-fluid 
model used here to describe superfluid neutron stars is 
practically equivalent to those use in previous studies, and 
we refer the reader to Andersson & Comer (2001a) and 
rrix et al. (2002| ) for a more extensive discussion about 
its physical motivations and justification. 

The plan of this paper is as follows: in Sect. ||we intro- 
duce the basic equations for the general two-fluid neutron 
star model, and in Sect. ||we develop its linear perturba- 
tion equations and show how to recover the special case of 
a single perfect fluid. In Sect. |4|we specialise to the simpler 
case of adiabatic oscillations of free, cold fluids, and we 
derive the necessary boundary conditions. In this section 



we also show that the separation into two distinct mode 
families is possible only for non-stratified models. Sect. |^ 
presents the numerical results concerning the background 
models, the eigenmode spectrum and its dependence on 
entrainment (resulting in avoided crossings), as well as the 
one-fiuid results where we recover the expected composi- 
tion g-modes. We present our conclusions and a discussion 
of necessary future work in Sect. |^. 

2. The general two-fluid neutron star model 

2.1. The general two-fluid equations 

The equations and notation for the Newtonian two-fiuid 
neutron star model used here are based on a more gen- 
eral formalism described in Prix (2002 ), which is the 



Newtonian a nalogue of a g enerally relativistic framework 
developed by [Carter (1989[) . In this section we will briefly 
summarise the general model and equations relevant for 
the present work, and we refer the reader to Prix (2002| ) for 
the derivation and more detailed discussion of this model. 

We describe a neutron star as a system consisting of 
two fiuids: a superfluid of neutrons, and a normal fluid of 
comoving constituents, which include protons, electrons 
and entropy (and generally further particles like muons 
etc). We denote the particle number densities for neutrons, 
protons and electrons as n^, np and ric respectively, and 
we use s for the entropy density. The velocities of the two 
fiuids are for the neutron fiuid, and for the fiuid 
of comoving constituents, the relative velocity A between 
the two fluids is therefore 



A = Vr 



(1) 



On the local ( "mesoscopic" ) level, a superfluid is con- 
strained to be in a state of irrotational flow, and its angu- 
lar momentum will be accommodated by a lattice of "mi- 
croscopic" vortices. However, for the global description of 
a neutron star we are more interested in a "macroscopic" 
description of the superfluid, obtained by averaging over 
volume elements containing a large number of vortices, 
but which are small compared to the dimensions of the 
neutron star. In this macroscopic description, the super- 
fiuid has a continuous vorticity and nearly behaves like an 
ordinary fluid (apart from small anisotropics due to the 
vortex-tension, which we will neglect). The absence of (lo- 
cal) viscosity still allows the superfluid to move relative 
to the normal fluid, but the presence of the vortex lattice 
now allows for a direct force between the two fluids. In 
the case of this force being zero, the vortex lattice moves 
with the superfluid and we refer to this situation as free 
fiuids. If the vortices are "locked" to the normal fluid (e.g. 
as can happen in the crust), the mutual force will be non- 
zero but strictly non-dissipative. This is known as "vortex 
pinning" . Only in the intermediate cases, where a friction 
force causes the vortices to have a different velocity from 
both the normal fluid and the superfluid, energy is dissi- 
pated. In this case the mutual force is usually referred to 
as "mutual friction" . 



Reinhard Prix and Michel Rieutord: Adiabatic oscillations of non-rotating superfluid neutron stars 



3 



An essential simplification of the present two-fluid 
model is that we neglect all electrodynamic effects, as we 
assume the charge densities of protons and electrons to be 
strictly balanced, i.e. 



Therefore we can effectively eliminate the electrons from 
our description, as their density and velocity is entirely 
specified by the protons. 

We note that in "tra nsfusive" models (as first set up 
in Langlois et al. (1998 )), i.e. models which allow for (3- 
reactions (n # p + e + S') between the two fluids, the total 
mass in the reaction has to be conserved in a consistent 
Newtonian description. Therefore we set 

TTib = rrin = nip + ■ (3) 

The respective mass densities of the two fluids can now be 
written as 



Pn ^rnhn^. 



and 



(4) 



and the total mass density p is simply p = Pn + Pc- 

The (local) kinematics of the system is completely de- 
scribed (up to arbitrary rotations and boosts) in terms 
of Pn, Pc, s and A^. The dynamics is determined by the 
internal energy density function £ or equation of state, 
which is a function of the form £ = ^(s, pn, pc, A^). This 
energy function defines the first law of thermodynamics 
for this system by its total differential, namely 



dS = T ds -\- p^ dpc + p" c?Pn + ce dA^ 



(5) 



This differential defines the dynamic quantities, namely 
the temperature T, the specific chemical potentials pi'^, pi", 
and the entrainment a as the thermodynamic conjugates 
of the kinematic quantities s, pc, Pn and A. The specific 
chemical potentials pt are related to the more common 
definition of the chemical potentials p via pt = p/nih- We 
note that p'^ is not simply the proton chemical potential, 
because adding a proton in this model implies adding an 
electron as well (due to (||)), and therefore one can see 
that p'^ = p^ -\- p^, where pP and p*^ are the respective 
proton and electron chemical potentials. 

The function a defined in (^) reflects the dependence of 
the internal energy on the relative velocity A between the 
two fluids, which characterises the so-called entrainment 
effect. This entrainment function a has dimensions of a 
mass density, and it will be useful in the following to define 
the two dimensionless entrainment functions £„ and Sc as 
_ 2a Pc 

= — = — £c ■ (6) 
Pc Pn Pn 

The "generalised" pressure P is introduced as the 
Legendre-conjugate of the energy density, namely by the 
usual thermodynamic relation 



— , and 



£ + P = pnp'' 



PcP 



sT, 



(7) 



which results in the total differential of the pressure func- 
tion P{p'',p'^,T, A2): 

dP = Pn dp" -I- Pn dp'' + sdT-a dA^ . (8) 



This generalised pressure P can be seen to reduce to the 
usual definition of the pressure of a perfect fluid in the 
case of A = 0. 

The equations of motion of this two-fluid system are 



derived from a "convective" variational principle in Prix 



(2002), and here we only present the resulting equations. 



The conservation of energy results in the following equa- 
tion: 

= /• A + /3rn, (9) 

where Pg and Pn are the creation rates of entropy and 
neutrons respectively, i.e. 

P, = dts + V ■ (sVc) , and P,, = dtPn + V • (pnVn) , (10) 

while the proton creation rate Pc has to satisfy Pc ~ —Pn 
for baryon conservation. The quantity /? in (|^) charac- 
terises the deviation from chemical equilibrium and its 
explicit expression if found as 



/3 ^ p^ - p" - - (1 



2£n) A^ 



(11) 



This is the Newtonian analogue of a result first found in 



the relativistic transfusive model by Langlois ct al. (1998). 
We note that there is an additional kinetic term with re- 
spect to the naive p'^ — p" in the case of relative motion. 
The mutual force density / (the sign convention is such 
that this force acts on the neutron fluid) is a direct in- 
teraction force between the two fluids, which in the case 
of a superfluid stems from vortex interactions like pinning 
or mutual friction. In order to ensure explicitly that the 
second law of thermodynamics, i.e. Pg > 0, is satisfied by 
(^, we can write the neutron creation rate Pn and the 
mutual force / in the form 

Tn = with S > 0, (12) 

f = 'qA + KxA, with 77 >0, (13) 

where the non-negative functions S and i] govern the 
beta-reaction rate and the friction force, while the vec- 
tor K allows for a non-dissipative Magnus-type force (i.e. 
orthogonal to the relative motion). A non-transfusive 
model, i.e. one that does not allow for beta reactions 
n ^ c, has S = 0, and free vortices correspond to 77 = 
and K = 0, while pinned vortices correspond to 77 = and 

The momentum equation for the superfluid neutrons 
is given by 

(at+-(;„-V)(t;n+enA)+V(p"+$)+e„AaV< = — /, (14) 

Pn 

while the equation for the normal fiuid reads as 

(dt + Vc ■ V)iVc - ecA) + V(p= + $) - £cAaV< 

+ -Vr=-l/ + (l-ec-e„)— A. (15) 

Pc Pc Pc 

The gravitational potential $ is related to the mass den- 
sities Pn and Pc via the Poisson equation 

^ AnGipn + Pc) . (16) 
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2.2. The static equilibrium background 

We consider a static background star, so we set 

Vn = Vp = A = 0, and r„ = Te = A = , (17) 

which by (|l^) implies the vanishing of the mutual force, 
i.e. / = 0. Because chemical reactions would not be neg- 
ligible on long timescales, we also assume the background 
star to be in chemical equilibrium, i.e. 



(18) 



These equilibrium conditions reduce the equation of mo- 
tion ([l4|) to 



(19) 



and with (|T^ this also implies that the background star 
is in thermal equilibrium, i.e. VT = 0. 

The static background has to be spherically symmet- 
ric, and therefore ([l^) and (^6|) lead to the following equa- 
tion for the background: 



fl"{r) + -Jl'ir) = -4^Gp(r) , 
r 



(20) 



where we have introduced the equilibrium chemical po- 
tential Jl = Jl'^ = Jl'^ , and the prime (') denotes the radial 
derivative d/dr. 

The equation of state allows one to relate the equilib- 
rium chemical potential JI directly to the total mass den- 
sity p at constant temperature, and therefore the back- 
ground is fully determined. The numerical method for 
solving this equation will be discussed in sect. ^. 

In the following it will be convenient to use the radius 
R and central density po of the static background as basis 
units for length and mass density, so the corresponding 
"natural unit" for frequencies is y/AnGpQ. All equations in 
the following are expressed in these natural units except 
otherwise stated. 

2.3. Entrainment and effective masses 

For small relative velocities A, we can separate the "bulk" 
equation of state from the entrainment by expanding 
f (s, pn, Pc', A^) in terms of A, i.e. by writing 



£ = (s, pn, Pc) + a^"^ [s, p„, Pc) A2 + OiA^) 



(21) 



where the quantities £^°^ and a^*^-* are evaluated at zero 
relative velocity A. The background equation of state f 
is therefore decoupled from the entrainment function a^"-* , 
and we can specify these two functions independently. The 
link between the entrainment function a and th e equiva- 
lent description in terms of effective masses m* (Andreev 
fc Bashkin 1975 ) has been discussed in previous work 
( Prix et al. 2002 ), and it can be shown that one can ex- 
press a in terms of the proton effective mass m* (which is 
generally a function of the densities), in the form 



2a 



p 

mb 



(22) 



The dimensionless entrainment functions^ £„ and Sc can 
then be expressed according to (ph as 



p 



and 



Pn 



(23) 



3. Linearised perturbation equations 

3.1. Oscillations of superfluid neutron stars 

We consider small perturbations with respect to the static 
equilibrium background described in Sect. 2.2. Linearising 
the equations of motion (^4|) and ( [l5| ) yields 

PndtiSv^ + SA) + p„V(5/I" + 5$) = Sf , (24) 
pMSvc - Ec SA) + pcV{SJl' + (5$) = -Sf - sVST, (25) 

where SQ denotes the Eulerian perturbation of the quan- 
tity Q. The perturbation of the mutual force ( |l3|) is given 

by 



Sf = rjSA + KX dA . 



(26) 



and the linearised energy conservation (^ and (^2|) to- 
gether with the condition of baryon conservation lead to 

dtSs + V ■ {sSv^) = 0, (27) 
at(5p„ + V-(pn<5-y„) = ^{5p'' -Sp!^), (28) 
dt Sp, + V- (pc Sv,) - -S {Sp' - Sp'') . (29) 

The perturbed Poisson equation ( p^ ) reads (in natural 
units) as 



V^(5$ = Spa + Spc 



(30) 



The system is closed by specification of the mutual force 
functions rj and k, the transfusion function S, and an 
equation of state which allows to express the dynamical 
quantities ST, Sp^ and Sp'^ in terms of the kinematic vari- 
ables Ss, Spn and Spc, thereby reducing the number of un- 
known perturbation quantities to 13, which corresponds 
exactly to the number of equations. 

3.2. The special case of normal-fluid neutron stars 

It is interesting to compare the superfluid neutron star 
case with the normal fluid case, where the two constituents 
n and c are moving together and form a single perfect 
fluid. This case is obviously just a subclass of the two-fluid 
case discussed so far, namely subject to the additional 
constraint = Vc, and thereforen 



Svn = Svc = Sv . 



SA = 0. 



(31) 



We note that iLindblom fc Mendell (200C| ) and more re- 



cently Andersson et al. (2002 ) have used a slightly different 
dimensionless function e to characterise the entrainment effect. 
The relation between e and Sc is given by e = £cPc/ (pn — £cp). 



Lindblom fc MendeU (1994, 1995) have imposed 5/3 = to 



recover the perfect fluid case. However, adiabatic oscillations 
of a perfect fluid only satisfy this condition in non-stratified 
stars (cf. sect. ~>.4), the constraint 5/3 — is therefore generally 
not met. 



Reinhard Prix and Michel Rieutord: Adiabatic oscillations of non-rotating superfluid neutron stars 



5 



By linking the two constituents together, the degrees of 
freedom have been reduced by three, and instead of the 
individual momentum equations (^) and (^5|), now only 
the sum of momenta can be required to be conserved, i.e. 



We point out that this system of equations is identical 
to the one used in Andersson & Comer (2001a), which 



was obtained from a Newtonian limit of the relativistic 



p [dtS-fi I Vi5$) I pryS'jT I pcVJ/I'' I jVJT-O. (32) 
We introduce the notation 

Pc = , p„ = x^p , and s^sp, (33) 

for the proton and neutron fractions Xc and x-^ and the 
specific entropy s, which allows us to rewrite the one-fluid 
equation of motion ( ^2|) in the slightly more familiar form 

'6P 
P 

where we used the fact that the total pressure differential 
dP of (||) in this perfect fluid case reduces to 

dP = p„ d/I" + p„ dp."" +sdT. (35) 

We can now compare (^2|) to the standard expression for 
the E uler equat i on of stellar oscil lations of non-rotating 
stars ( Cox 198C ; Unno et al. 198£ ) , which is usually writ- 
ten as 



equations. It is also related to the equations of Lindblom 
& Mendell (1994), which are expressed in the "orthodox" 



dtSv + V 



5$ 



[6(3 Vxc + (5r VS] = , 



(34) 



dtSv 



p 



-C(r) (V ■Sv)A = 0, 



(36) 



where C{r) is a function of the background, and A is the 
so-called Schwarzschild discriminant that is responsible 
for presence of g-modes. By comparing eqs. ( ^2| ) and (^^, 
we see that A will be non-zero (indicating the presence 
of g-modes) whenever Vs^ 7^ or Vxc 7^ 0. This reflects 
the well-known fact that any type of stratification, either 
in specific entropy J or in the chemical composition x. 



will result in g-modes, as pointed out by Reisenegger & 
Goldre |ich (199^ ). 



4. Adiabatic oscillations of "free", cold fluids 

4.1. Reduction to a ID eigenvalue problem 

In order to close the system of perturbation equations 
(p^)-(pO|), we need specific models for the mutual force 
/, the transfusion S, in addition to an equation of state 
of the form £{pn,Pc, s), all of which are highly dependent 
on microphysical models and are rather poorly known at 
the present stage. For this reason, we will postpone the 
inclusion of these effects to future work, and focus on the 
case of purely adiabatic oscillations (i.e. S = 0) of free 
fluids (meaning / = 0). We will further neglect tempera- 
ture effects (which is generally a very good approximation 
except for very young neutron stars), so we set s = and 
r = 0. The resulting simplified system of equations is 



formulation of superfluids (Landau & Lifshitz 1959), while 
the present description is based on the "canonical" ap- 
proach introduced by [Carter (1989 ) . 

Using the equation of state we can link the density 
perturbations Spx to 6JP^ (with the constituent index no- 
tation X = n, c) to linear order, namely 



Spn Snn (5^" + '^nc Sp" 
6pc = Sen (5/1" + Sec Sp" 



(42) 



where the symmetric "structure matrix" Sxy is defined as 



'XY 



X,Ye{n,c}. (43) 



dpY dpx dpY ' 

Due to the dual role of the pressure P in (|^) and (^, we 
can equivalently express Sxy as 



Sxy — 



dpx 



d^P 



(44) 



dtdpn + V ■ jpnSVn) = 0, 



dp^ dp^ dp^ ' 

We note that although we have assumed free fiuids, i.e. 
there is no direct force acting between them, the fluids 
are nevertheless locally coupled by the equation of state; 
we can distinguish two sources of this coupling, one is due 
to the non-diagonal term S^c in (^) , while the second is 
due to the entrainment terms £x- 

The background quantities Jp^ can be seen in (^) to 
behave like the gravitational potential <I>; this means in 
particular that their gradient is always finite, even at the 
surface. Therefore SJP^ is finite everywhere, contrary to 
6px which can diverge at the surface when p' —^ —oo. 
This is seen from the relation Ap — Sp + between 
the Lagrangian perturbation Ap and the Eulerian Sp, for 
a radial displacement ^'^ . On physical grounds Ap must be 
bounded everywhere (as it reflects the physical property 
of a fluid element), while the flrst-order Eulerian quan- 
tity dp diverges at the surface whenever p' — > — oo and 
7^ at r = i?. This might seem problematic for the va- 
lidity of the equations, but it only reflects the fact that in 
this case even an infinitesimal displacement of the surface 
will lead to a finite (as opposed to infinitesimal) Eulerian 
density change there. By considering Lagrangian instead 
of Eulerian variables, it can be shown that the physical 
solution is still well behaved even if — > oo at the sur- 
face. In this case the first-order quantity Sp no longer ap- 
proximates the physical Eulerian density change, but the 
divergence is such that the Lagrangian first-order quan- 
tity Ap is still perfectly regular. If one wanted to impose 
that dp should be bounded everywhere (as did Lindblom 
(37) fc MendeU (199^)), then this situation would be inverted 



dtSpc + V • (pc Svc) 

dt [(1 - £n) + SnSVc] 

dt [(1 - Ec) + ScSvn] 



0, (38) 

-V((5p"-l-(5$) , (39) 

-V {Sp'' + 5<^) , (40) 

Spn + Spc, ■ (41) 



and the Lagrangian quantity Ap would diverge, which is 
unphysical indeed. 

From a numerical point of view it seemed better to 
solve directly for the well-behaved Sp?^ instead of the po- 
tentially diverging Spx, by using (B^) to substitute for Spx- 
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We note that the coefBcicnts Sxy in this expression will 
generally diverge (or vanish) at the surface, depending on 
the equation of state and reflecting the behaviour of Spx- 
The system of equations (p^Bl]) for eigenmode solutions 



(45) 
(46) 
(47) 
(48) 
(49) 



where ©(r') means of order r' or higher. Another require- 
ment is the regularity of the solution at singular points of 
the equations. From eqs. (pZ]) and (pS) we see that Wx 



of the form dQ{x, t) = SQ{x) e*'^* now yields 
V((5aI'= + (5$) = -iuj[{l- ec)Svc+enSvn] , 



must vanish as 0{r) at the origin, while from eqs. (59) 
and (|^) we can derive the further regularity requirement 



(63) 

requires 
(64) 



where we have introduced the convenient "structure vec- 
tor" kx, which is defined as 



E 

Y=n.c 



XY 



X e {n, c} . 



(50) 



For a spherically symmetric background, we can separate 
the radial and angular dependence and obtain solutions 
with definite quantum numbers I and \m\ < I using the 
ansatz 



6^{r,d,ip) = S'^{r)Yi"'{9,ip), 
SJP'ir^e,^) = SlP'{r)Yr{9,^), 



(51) 
(52) 



^1 + Wx + rW^ - l{l + l)Vx ^ 0{r^) . 

Regularity of solutions of Poisson's equation (|6 

5^-^0{r^) for l>0, 

while for radial oscillations {I = 0) only (5$' ~ 0{r) is 
required. These constraints are automatically satisfied by 
(p2) for I > 2, but they are stronger requirements than 
(p2) in the cases I = 1 and I — 0. 

4.2.2. At the surface 

At the outer surface (r = R) we need to ensure the conti- 
nuity of the gravitational potential $, which results (e.g. 
see Ledoux & Walraven (195S)) in the boundary condition 



R 



-<5<f (i?) = -exi?-) 



(65) 



Vx{r)S^-^Ux{r)T^, (53) 



where the spherical harmonics ¥"'■{6, ip) are the eigenfunc- 
tions of r^V^i^™ = + 1) F;", and fi, S and T form 



the orthogonal harmonic basis, defined as 



rrvr, 



grn ^ 



TV' = V X il , 



(54) 



see Rieutord (1987 ) for details. The three-dimensional 
eigenvalue problem (|4^)-(|49|) has now been reduced to the 
following one-dimensional problem: 



(rpnW^n)' 


-;(Z + l)p„K 




+ SncSJl'] , (55) 


(rpcWc)' 


- l{l + l)pcl4 


= —iur'^ [iSnc'^M" " 


f 5cc5m1 , (56) 




r {SJi'" + S^') 


= iu; [(1 - £„)iy„ 


+ EnW,)] , (57) 




r {6Tf' + <5$') 




fecWn)], (58) 




^5p^ - 5^ 




- EnV,)] , (59) 






= iuj[{l~ec)Vc + 


ScVn)] , (60) 


{rH<^>') 


' -l{l + l)d<i> 


= r^iknSJl^ + kcSJl"). (61) 



where is the radial displacement of the surface, 
and p(R—) is the inner limit of p{r), i.e. p{R—) — 
limr_>i{_ p{r) . In the present work we will only consider 
stars with vanishing density at the surface^ i.e. p{r) 0. 
The conservation equations ( p5[ ) and ( |5^ ) contain a (regu- 
lar) singularity at the surface r = R, because p' / p diverges 
as p — !■ 0. Therefore we need an additional regularity con- 
dition that can be obtained by a Frobenius-type expan- 
sion around the surfacc^,i.e. /i = 0. We first rewrite these 
equations as 

Pxw , ■ 2T.SxY5'il^^ 
'—^Wx + it^T — 



{rWx)'-l{l+l)Vx = 



Px 



(66) 



and using the definition ( p3| ) of the matrix Sxy together 
with the equilibrium condition pi"' = = /I', we can 
write 



P'x 



kxp! : 



(67) 



The axial velocity component Ux is decoupled and cor- 
responds to zero frequency w = 0, therefore all non-zero 
frequency eigenmodes are purely polar. We note that the 
horizontal velocity equations (^9|) and ( |60| ) only hold for 
I > 0, because Sq = 0. 

4.2. Boundary and regularity conditions 
4.2.1. At the center 



where kx has been defined in (|5C|). In order to simplify 
the analysis, we will assume that the equation of state 
becomes decoupled in the limit of vanishing density, i.e. 



P(a2",a2^)^P„(/I") + Pc(m^): 



(68) 



such that pn ~ /5n(/Lt") and Pc '■^ pdp'^), and kx ~ Sxx in 
this limit. In this case we can write the asymptotic form 
of M) as 



S 



{rWxY - Kl + l)^x = -— [rpWx + iuor^Sp 



(69) 



It can be shown that the representation (plD-([53D of a reg- 
ular physical quantity requires the following asymptotic 
behaviour of the radial functions as r — > 0, 



Sp^ S<S> Wx ^ Vx ^ 0{r^) , Xe{n, c}^ 



(62) 



We note that [Lindblom fc Mendell (1994 |1995[ ) have set 
the right-hand side of (|65|) to zero in their homogeneous model, 
which is inconsistent if one allows for surface displacements . 

* We can conveniently use p as an effective radial coordinate 
at the surface, because it is well-behaved and monotonic (i.e. 
p,{R) = and p'iR) < 0). 
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where the expression in brackets is regular for Jl —^ 0. 
The behaviour of Sxx/px in this limit can be analysed 
by writing px in the asymptotic form px ~ (/J"^)" (with 
positive a), which is possible due to the assumed decou- 
pling. Therefore we find Sxx = dpx/djj?^ ~ 
we have identified the singularity in (|6^) as 

^xx ~_i 



/i" ^, and so 



PX 



(70) 



Regularity of the solution therefore requires the following 
asymptotic behaviour at the surface: 

r'jl'Wx + iiur^SJl^ ~ 0(jl) , for X e {n, c} . (71) 



4.3. Decoupling "ordinary" and "superfluid" modes? 

In this section we discuss a change of variables that has 
been used in several previous studies of oscillations of su- 



perfiuid neutron stars ( 


Lindblom & MendcU 1994, 


1995, 


200C ; ISedrakian & Wasserman 2000; 


Andcrsson & Comer 


2001a), namely 



SA = SVc — SVn , 
Sv = Xc SVc + Xn Sv^ . 



5(3 = - Sp" , 

Sp = Xc Sp"^ + Xn SJl'^ 



(72) 



This choice of variables is motivated by the intuitive idea 
that the additional degrees of freedom of a second fluid 
should allow for two different types of motion, charac- 
terised roughly by the two fluids being either "co-moving" 
or "counter-moving" , and which are sometimes referred to 
as "ordinary" and "superfluid" modes. By choosing such 
"adapted" coordinates ([z^), one might hope to separate, 
or at least simplify the system of equations, but we will 
see that this is generally not the case. 

Using the relations (^ ) and the definition ( |67|) of /c„ 
and kc, and defining k = + kc, we can express the 
density perturbations Spx as 



5pn = k^Sp 4 
5pc = kc5p + 

and therefore 



(5„ 



-5„ 



Xckn) 5/3 , 
+ Xnkc) (5/3, 



Sp = kSp — (Xckn — Xnkc) 5(3 . 



(73) 



(74) 



In the case of a spherically symmetric background consid- 
ered here, we can use (|67|) to express 



Xnk^ 



-k— X 
P 



(75) 



and using these variables and relations, we now rewrite 
the eigenmode equations (^^-(^Tj) in the form 



kdt5p + V ■ (pSv) 
dtSv + V{Sp, + S^) 
y'^6^-kSp 

ICdtSp + V ■ (pXnXcSA) 
2ap' 



k^dt5(3 
. P 



x'^ [5(3 e"- 



k^5(3 
. P' 

k — dtSp + pSv^ 
L P 



1 



PnPc 



dtSA + V6(3 



0, 



(76) 
(77) 
(78) 

(79) 

(80) 



where e'' is the radial basis vector, and we introduced the 
abbreviation /C = x^kn + x^kc — Snc- 

We see that the "ordinary" -type of motion {Sp, Sv, 
S^) does not decouple from the "superfiuid"-type vari- 
ables {SA, S(3) whenever there is stratification, i.e. when 

7^ ! This can be understood as follows: while a non- 
zero relative velocity SA can be regarded as a characteris- 
tic of a superfiuid mode (as opposed to modes in a single 
fluid), the chemical equilibrium deviation S(3 is generally 
non-zero even for a single (but non-barotropic) fluid. In 
a stratified fiuid, any general adiabatic motion will drive 
a fluid clement out of equilibrium, i.e. S(3 nonzero is not 
characteristic for eith er "superfluid" or "ordinary " mod es 
(contrary to claims in [Lindblom fc MendcU (1994 |l995|) ), 
it is a general feature of modes in stratifled fluids, and 
therefore the choice of variables (|7|) does not lead to a 
decoupling of the system in this case. 

However, it is interesting to consider for a moment this 
special case of a non-stratified background (which proba- 
bly never applies in real neutron stars) . Setting the proton 
fraction Xc to a constant, we can separate the equations 
into two decoupled sets. One system describes "ordinary" 
modes, namely 



kdtSp + V • {pSv) = , 
dtSv + V{Sp + S^) = 0, 
V2(5$ -kSp = 0, 



(81) 
(82) 
(83) 



which are seen to be independent of the entrainment e 
as well as of the coupling through the "bulk" equation of 
state, i.e. iSnc- The second system of equations governs the 
"superfluid" modes and reads as 



k 



Sr. 



XyiXq 



1 



dtS(3 + V -{pSA) 
2ap 



0. 



PnPc 



dt SA + VS(3 = , 



(84) 
(85) 



We see that contrary to the ordinary modes, the superfluid 
modes do depend on the coupling through entrainment a 
and the equation of state, i.e. S^c, but they are completely 
decoupled from the gravitational perturbation 5$, as they 
leave the total density unchanged, i.e. Sp = 0. 

We have therefore shown that in the non-stratified 
case there exist two separate families, namely "ordi- 
nary" modes {Sv, Sp, S^, 0, 0) and "superfiuid" modes 
(0, 0, 0, (5A, (5/3). One of our numerical models (see next 
section) has a constant proton fraction a;c, and we will see 
the present analysis confirmed by the numerical results for 
this model. In the general stratified case, however, these 
two mode families are coupled and such a clearcut sepa- 
ration is not possible. 

5. Numerical results 

5.1. Equation of state: two-constituent polytropes 

We use a simple class of two-constituent equations of state 
which is very convenient to explore the properties of a 
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two-fluid system, namely the following "generalised poly- 
trope", defined as[] 



Pc 



Kn pI" + Kc p^" 



(86) 



which simply consists of the sum of two ordinary poly- 
tropes. For regularity of the chemical potentials (||) in the 
limit Pn and pc — > 0, the polytropic indices must sat- 
isfy 7n > 1 and 7c > 1- This equation of state allows the 
explicit inversion 



pnijin 



Kn7n 



and Pcip-"") 



P" 



Kc7c 



(87) 



where we have introduced iVx = l/(7x ~ !)• We see that 
in chemical equilibrium, i.e. Jl^ — Jl'^, the two fluids share 
a common outer surface. The equilibrium proton fraction 
Xc can be expressed as 



Xc = { I 



(7c^c)^- 
(7nKn)^» 



■P 



which shows that the proton fraction is constant whenever 
7,1 — 7c , while the behaviour in the case of different indices 
falls into the two categories: 



< Xr > 1 



1 < Xc — > , 



for 7„ < 7c , 
for 7n > 7c ■ 



(89) 
(90) 



5.2. Calculating the background models 



The equilibrium background solution is determined by 
equation (pO|), together with the regularity requirement 
/i'(0) — 0, and the boundary condition of vanishing pres- 
sure at the surface, i.ej^ /I(l) = 0, where the surface of 
the static background star is situ ated at i? = 1 in the 
natural units defined in Sect. 2_^. Using the equation of 
state we can express p = p(jl) in chemical equilibrium. 



and therefore equation (20) can be written as the follow- 
ing nonlinear eigenvalue problem. 



P" + -P = ^P{p) , 
r\ 



done by solving in step fc+1 the following linear eigenvalue 
problem, 



pU 



~l^'k+l - 



fe+l — Pk+l 

\ Pk J 



(93) 



where pk is the solution of the of the previous step 
k. The equation in each step is solved using the spec- 
tral linear eigenvalue solver package LSB developed by 
L. Valdettaro and M. Rieutord. With a resolution of 40 
Chebychev polynomials and a Gauss-Lobatto collocation 
method, this iteration converges to about machine preci- 
sion (i.e. ~ 10~^^ relative difference between successive 
steps) in about 20 steps. Another practical advantage of 
this method is that we calculate the background quantities 
on the same Gauss-Lobatto grid we use for the numeri- 
cal code for the eigenmodes. Even when using a different 
resolution for the eigenmode-calculation, the Chebychev 
expansion provides a canonical C°° interpolation which al- 
lows us to easily "re-grid" the background solution. 

In principle we can solve the background for any given 
equation of state, but we restrict our attention in this work 
to the class of two-constituent polytropes (|8^).The reason 
for this choice is their analytic simplicity and because our 
main emphasis is to clarify the qualitative properties of 
superfluid neutron stars rather than to construct a physi- 
cally very precise model. This would in any case be quite 
impossible in a Newtonian description because of the ne- 
glect of relativistic effects, in addition to the important 
uncertainties in our current knowledge of the equation of 
state of neutron stars. 

In the numerical analysis we consider the three differ- 
ent background models deflned in table |^, and which are 
represented in Fig. |l|. These three models correspond to 
three different types of behaviour at the surface. In the 
case of Model I one can easily flnd the backgro und so- 
lution analytically, as shown in Prix et al. (2002 ), which 
allows us to check the numerical method of calculating 
the background, and we find a maximal relative error of 
10~^^ between the numerical and the analytic solution for 
model I. Model II represents a generic "stiff" model sim- 



(91) ilar to those used in Comer ct al. (1999 ) and Andersson 



where the (dimcnsionlcss) eigenvalue A is given by 
P[Po) 



(92) 



The eigenvalue A represents the a-priori unknown radius 
of the star, and thereby the actual values of our natu- 
ral units for a given central density pQ. The method used 
here to solve this equation is to iterate a linear eigenvalue 
problem that converges to the solution of (|9l|). This can be 



This equ ation of state has been used p reviously by Comei 
et al. (1999) and Andersson et al. (2002) to study two-fluid 
oscillations in general relativity. 

^ Strictly speaking p is determined only up to a constant, 
which i s usually fixed such that vanishes together with the 
pressur i. 



et al. (2002 ), which has infinite density gradients at the 
surface. Model III is of a "soft" type with vanishing den- 
sity gradients at the surface. These different types of be- 
haviour at the surface are quite analogous to the case of 
the usual one-constituent polytropes for different poly- 
tropic indices. We note that model I is the only non- 
stratified model, because it has a constant proton frac- 
tion Xc (and finite density gradients at the surface) , while 
models II and III have a non-zero composition gradient, 
i.e. x'c ^ 0, as expected from (|8^). 

For easier comparison of the frequencies given in the 
next section in units of y/inGpo, we provide in table |^ 
the conversion factors into three important systems of 
units, namely the SI unit Hz, the " Cox" units \/GM/R^ 
(variants thereof, like those used by Lindblom fc Mendel] 



(1994 ), only differ by a constant factor) and the "geo- 



0.2 



0.15 



0.05 
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Model I Model II Model III 

' I ' I ' 



0.2 



0.15 




0.05 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 
r/R r/R r/R 

Fig. 1. Density profiles and proton fraction Xc of the background models I, II and III as defined in table |l|. po is the 
central density. In this figure wc have divided the neutron density by 5 in order to obtain similar magnitudes for the 
different curves. 





7n 


7c 


K " 




PO 
Pnuc 


a;c(0) 


AI 

Mq 


_R_ 

km 


I 


2.0 


2.0 


0.01 


0.09 


10.83 


0.10 


1.414 


10.7 


II 


2.5 


2.1 


0.01 


0.20 


3.20 


0.10 


1.440 


14.3 


III 


1.9 


1.7 


0.01 


0.09 


18.65 


0.10 


1.412 


9.4 



Table 1. Parameters describing the background mod- 
els I, II and III, based on the two-constituent poly- 
trope (^). Further shown are the central mass den- 
sity Po, the central proton fraction Xc(0), the total 
mass M and the radius R. ° The units of the coeffi- 



cients Kx are c^pnuc^' 



'''^j where pn 



i.eexio^-^gcm-a 





I 


II 


III 


kHz 


38.8318047 
3.14159265 
0.270521502 


21.1116405 
2.61829313 
0.149723677 


50.9444446 
3.40801371 
0.354013053 


e'IGM 



Table 2. Conversion factors for our freq uency units 
(\/47rGpo) hito different systems of units (SI, |Cox (1976 ) 
and "geometrised" units) for models I, II and III. 



metric" units (? jG M typically used in general relativity 
( [Comer et al. 1999| ; [Andersson et al. 2002] ). 



5.3. The two-fluid oscillation modes 

The eigenmode equations ( p7|)-(|6T| ) together with the 
boundary conditions of Sect. |4.^ form a linear eigen- 
value system which we solve numerically using the spectral 
solver of the LSB-package. The convergence of the results 
was determined by increasing the resolution starting from 
40 Chebychev polynomials up to 80, and we found the 
changes in frequency decrease very quickly to about 10~^ 
(or better), which is why we give the frequencies with 
nine decimals in tables || to ||, corresponding roughly to 
the convergence achieved by the numerical method. This 



frequencies represent a physically measurable prediction 
in any sense. 

5.3.1. Eigenmodes of "locally uncoupled" fluids 

In this section we consider the case of zero entrainment, 
i.e. a = and £n = Sc ~ 0. We refer to this situation 
as "locally uncoupled" fluids, as it is important to note 
that the two fluids are nevertheless coupled "globally" 
through the perturbation of the gravitational potential 
and (|6l|). We consider the cases of radial (Z = 0), dipolar 
{I = 1) and quadrupolar {I = 2) oscillations, which differ 
qualitatively in some properties and boundary conditions 
(see Sect. 4.2), while all higher I cases are qualitatively 
very similar to Z — 2. The lowest eigen- frequencies for 
these three values of I are shown in tables H, ^ and || respec- 
tively. We label these modes in analogy to the one-fluid 
case as f- and p- modes, and group them in pairs where 
the lower frequency mode is labelled as "o" and the higher 
frequency one as "s" . The pairs of p-modes are indexed in 
the order of increasing frequency. We emphasise that this 
labelling is a pure convention, as one can generally not say 
that these modes would be either co- or counter- moving, 
or that the subscript would exactly reflect the number of 
radial nodes. 

Let us flrst consider the special separable case of the 
non-stratified model 1. The first three pairs of eigenfunc- 
tions are presented in Fig. ^ and we see that in the "o" 
modes the two fluids are comoving, resulting in a non-zero 

and they also remain in strict chemical equilibrium, 
i.e. dp = 0. These "ordinary" modes are actually iden- 
tical to the normal-fluid modes of the same background 
(see Sect. 5^). In the case of the s-type modes the two 
fluids are counter-moving in exactly such a way that the 
total density remains constant, i.e. 6p — and therefore 

= 0, while the fluids are driven out of chemical equi- 
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0.2 0.4 0.6 0.8 
r/R 



0.4 0.6 0.8 

r/R 



0.2 0.4 0.6 0.8 

r/R 



Fig. 2. The first three quadrupolar (/ = 2) eigenmode doublets (o)/(s) for the background model I. In the case of 
the o - modes, the two fluids are comoving to numerical precision, therefore the 5v,^ and 6vc curves can not be seen 
separately. 



librium, i.e. 5j3 7^ 0. The number of radial nodes in 5v^ is 
the same for the o and s modes, and corresponds exactly 
to their index. All these results confirm the analytic pre- 
dictions for non-stratified models in Sect. |4.3| . However, it 
would be wrong to assume that these properties are gen- 
erally true for superfiuid oscillations. Stratification makes 
this picture more complex, even in the case of locally un- 
coupled fluids considered in this section. If we look at the 
first three pairs of eigenfunctions for the model II in Fig. ^, 
we see that the "o" modes are not comoving at all (only 
the P is nearly comoving), and they have non-zero (5/3, 
while the "s" modes have non-zero 5p and (5<f>. One can 
not say either that the relative amplitude of 5(3 would 
be different between the o- and s- cases, as Lindblom & 
Mende [r(199^ ) wrongly induced from the properties of the 
f°-mode, the only eigenmode they presented. 

In the case of the low-order modes presented in Fig. H, 
the number of radial nodes still seems to correspond to 
the index, and the fluids are roughly in opposite phase in 
the s- modes, while they are approximately in phase for 
the o-modes. Even this, however, is not true in general, 
as can be seen in Fig. ^ which shows the higher order pg- 
modes. In this case neither o- nor s- are dominantly in 
phase or in opposite phase, and the two radial velocities 
have different numbers of radial nodes. Therefore neither 



I = 




UJII 


i^iii 


f° 


0.616 801012 


0.860 501 159 


0.539 820 916 


f 


0.825 395 141 


1.004 218 360 


0.713 202 951 


p? 


1.272153 763 


1.650 440 676 


1.114 032 342 


v\ 


1.398 557067 


1.780 808 966 


1.186 350 906 


P2 


1.855 852 617 


2.326 264 710 


1.582 750 279 


Pl 


1.949 822 942 


2.573 698 536 


1.675 075 045 


P3 


2.418 457 671 


2.985 429110 


2.018 290 454 


Pl 


2.493 326179 


3.352 297 316 


2.169 556 932 


P4 


2.970 977 248 


3.638 960 946 


2.445 611455 


Pl 


3.033169 591 


4.122179 788 


2.658 934 809 



the inc ex nor the o/s label bear any reliable information Comer (2001a ). 



Table 3. Frequency spectrum of radial eigenmodes (Z — 0) 
for models I, II and III in natural units wq = V^^rGpo- 



about the properties of the modes. This behaviour is pos- 
sible because the that this eigenvalue problem is not of 
Sturm-Liouville type except in the non-stratified case. 
An interesting fact to notice in tables || to || is that the 
fundamental modes (f° and f^) are the lowest frequency 
modes in the spectrum, in other words there are no g- 
modes present (which usually lie far below the f-mode) in 
these s uperfluid m odels. This confirms the n umerical find- 
ings by Lee (1995 ) and the local analysis of Andersson fc 
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0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 1 
r/R r/R r/R 

Fig. 3. The first three quadrupolar (/ = 2) eigenmode doublets (o)/(s) for the background model II. The value of 5(3 
in the lower left figure has been multiplied by 100, in order to show that it is non-zero. 



I = 1 








I = 2 




UJII 


i^iii 


f° 


0" 


0" 


0" 


f° 


0.390 550 961 


0.424 294 338 


0.376 662 787 


p 


0.389 835 134 


0.440 176 989 


0.387 577 390 


f 


0.526 990 499 


0.604 904 572 


0.516 636 947 


p? 


0.898 011966 


1.191707 859 


0.798 353 611 


P? 


1.101827434 


1.423 800 575 


0.980 798 266 


p! 


1.040 570 747 


1.280 907 920 


0.895 985 737 


v\ 


1.206 881695 


1.514 568 856 


1.039 415 772 


P2 


1.518 621841 


1.930 140 941 


1.317 638 486 


P2 


1.723 444 868 


2.164156 378 


1.478 705 720 


Pl 


1.621023 028 


2.105 369 668 


1.384 086 979 


Pl 


1.806 873 473 


2.380 237199 


1.557 015 636 


P3 


2.099 218 641 


2.609 579 518 


1.770 672 472 


P3 


2.310 315 782 


2.858 007106 


1.932 213 886 


vl 


2.179 811073 


2.908 118 483 


1.887 282 252 


Pl 


2.379 342 400 


3.200 741221 


2.072 313 682 


P4 


2.662 623 840 


3.274 392 031 


2.206 080 421 


P4 


2.879 785 468 


3.533 905 837 


2.372 551 143 


Pl 


2.729 000 908 


3.692 629 614 


2.385 023 708 


Pl 


2.938 448 588 


3.997 594 591 


2.576 350 147 



Table 4. Frequency spectrum of dipolar eigenmodes (/ = 
1) for models I, II and III in natural units loq = ^/AirGpa. 
" The zero frequency modes correspond to the analytic re- 
sult of a constant displacement field ^. This solution can 
not be produced by our code because the equations have 
been expressed in terms of velocities instead of displace- 
ments. 



Table 5. Frequency spectrum of quadrupolar eigenmodes 
{I — 2) for models I, II and III in natural units = 
\/47rG>o- 



frequency (see [Dintrans fc Rieutord 200l| ; [Rieutord fc 



The absence of g-modes can be made clearer when 
acoustic modes and surface gravity modes are filtered 
out. The latter modes are easily removed by suppressing 
surface motions and imposing therefore Wx = at the 
star surface. Acoustic modes, on the other hand, are fil- 
tered out by using the so-called anelastic approximation 
which makes an expansion in powers of the Brunt- Vaisala 



Dintrans 2002 ). Using this approximation mass conser- 
vation now reads 

V • (p„(5t;„) = and V-{pc5v,)^{) (94) 

instead of (^ ) and (^). In the case we have considered, 
i.e. that of no entrainment (a = 0), the equations of mo- 
tions read: 



V(5/I" -t- (5$) = -iLoSvi^, 
V{5'jl'' + (5$) = -ibjSvc , 



(95) 
(96) 
(97) 
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e =-0.5 



e = -0.2 



e = 0.1 




Fig. 4. The pe eigenmode doublet (o)/(s) (for I ~ 2) for 
the background model II. The radial velocity Sv^^ has been 
multiplied by 10 for better visibility. Neither mode is co- 
or counter-moving, and the number of radial nodes differs 
between the two fluids. 



Using (p^, we eliminate the velocities and are left with 
the system: 



<5$)) = 0, 
S^)) = 0, 

fc„ SJl"" 



fee Sfl^ 



(98) 
(99) 
(100) 



where we see that the mode frequency has disappeared. 
Using the boundary conditions, it turns out that the only 
solution is SJl'^ = 5/1" = 6^ — 0, thus showing that no 
eigenmode exists when acoustic and surface gravity modes 
are filtered out. 

However, the presence of g-modes due to chemical 
composition gradients in normal-fluid neutron star models 
has been pointed out by Reisenegger fc Goldreich (1992| ), 
and their possibly observable excitation in a coalescing 



binary neutron star has been discussed by Reisenegger & 
Goldre |ich (l"994| ) and |Lai (1994] ). We wiU see in Sect. ^ 
that these predicted composition g-modes do indeed ap- 
pear in the normal-fluid case. In principle the presence 
or absipnce uf Lhese luudcs could Lherefutc be used as a 
possibly observable liidicaLur foi' supcrfluldiLy in neuLi'uii 
stars. 
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Fig. 6. Avoided crossing between the p^ and the p\ modes 
of Model II (for 1 — 2), indicated by the box in Fig. |[ 
The three columns show the corresponding eigenfunctions 
for three different values of entrainment. The small arrows 
indicate the "evolution" of the eigenmode when increasing 
e. The Sv^ curves in the two lower left figures have been 
magnified by a factor of 10. 



5.3.2. The effect of coupling by entrainment 

In this section we study the dependence of the mode- 
frequencies and properties on the coupling by entrain- 
ment. Obviously, we only need to specify one entrain- 
ment function, £c say, as En is then determined by (^. 
Because the uncertainties and differences of the "realistic" 
models for Sc provided by nuclear physics calculations so 
far are still considerable, we chose the simplest entrain- 
ment model, namely Sc = s being a constant. The value 
of this constant e can be related to the proton effective 



mass m* {Uih, and is roughly constrained from the nu- 



clear physics calculations (Chao et al. 1972; Sjoberg 1976 
p3aldo et al. 1992 ; Borumand et al. 1996) to lie in the range 



0.3 < e < 0.7. We nevertheless consider the broader range 
between —0.8 < e < 0.8 to demonstrate the qualitative 
behaviour more clearly. This will also show that the "lo- 
cally uncoupled" case e = (considered in the previous 
section) is not special, contrary to what one might have 
expected. 

The results for the mode-frequencies as functions of e 
for the three background models are represented in Fig. |[ 
In the case of the non-stratified model I, we observe the 



predicted (sect. 4.3) decoupling, and in particular the in- 
dependence of the "ordinary" -type modes of entrainment. 
Because of this decoupling the respective frequencies of 
the two mode families can simply cross each other when 
e is varied. In the generic stratified models (model II and 
III), the modes of the doublets are coupled and avoided 
crossings result when mode-frequencies come to close to 
each other, as also found recently by Andersson et al 



(2002 ) . In this process of avoided crossing the two modes 
seem to exchange some of their respective properties of 
being dominantly "co-" or "counter-moving" , as can be 
seen in Fig. P), and they also can exchange their number 
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Fig. 5. Eigenmode frequencies uj (for Z = 2) as functions of entrainnicnt e for models I, II and III. The full and dashed 
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Fig. 7. Avoided crossing between the pf and the P2 modes 
of Model II (/ = 2), indicated by the hatched box in Fig. ||. 
The three columns show the corresponding eigenfunctions 
for three different values of entrainment. The small arrows 
indicate the "evolution" of the eigenmode when increasing 

£. 



of radial nodes, as we see in the avoided crossing of the p^ 
and P2 in Fig. 0. 

Another important conclusion can be drawn from 
Fig. |5[ namely that the "locally uncoupled" case £ = 
discussed in the previous section does not represent a spe- 
cial case in any respect, because the two fluids are always 
coupled through (5$. The efl'ect of e is simply to change 
the coupling, but no configuration is completely uncou- 
pled. On can see in Fig. || that several avoided crossings 
happen practically at £ = 0, which is the case in particular 
for the p6-modes of model II presented in Fig. ^. 



I = 2 


UJl 


UJII 


UJlII 


g4 




0.012105 268 


0.011489 709 


g3 




0.015 003 335 


0.014 157 505 


g2 




0.019 814 997 


0.018 492 575 


gl 




0.029 631 110 


0.026 880 058 


/ 


0.390 550 961 


0.424 310 492 


0.376 717 911 


Pi 


1.101827434 


1.477 988 230 


0.988 324 062 


P2 


1.723 444 868 


2.348 478 094 


1.533 337 250 


P3 


2.310 315 782 


3.163 853 143 


2.050 348 013 


P4 


2.879 785 468 


3.954 289 860 


2.552 745 872 



Table 6. The eigenvalue spectrum (for / = 2) for the 
models I, II, and III. The non-stratified model I has no 
g-modes. 



5.4. The one-fluid case: recovering the g-modes 



Following the discussion in sect. 3.2, the one-fiuid case is 
defined by Sv^ — Svc = Sv. We only have one Euler equa- 
tion in this case, which in the harmonic decomposition 
( p4|) has the two components 



Xn SJl^ + Xc SJl'^ + (5$ 



6^' = -'-^W, 
r 



-icuV . 



(101) 
(102) 



These two equations replace (p7|)-(|60D, while the remain- 
ing equations of this system are unchanged (subject to 
the substitutions Wn = Wc = W and Vn = Vc = VY 
The eigen-frequencies of this system are shown in table || 
where we see the presence of composition g-modes as ex- 
pected for all models wit h stratification. This is consistent 
with the prediction by Reisenegger fc Goldreich (1992 ) 
and the numerical findings of Lee (1995| )^ 
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small radial velocity Sv^ at the surface, while it is very large in the case of p-modes. 



We also see that in the non-stratified model I, the 
one-fluid frequencies and modes correspond exactly to the 
corresponding "ordinary" -type solutions of the two-fluid 
case (see table |^ and Fig. |^), as would be expected from 
the separabiUty of the system as discussed in Sect. |4.3| . 

We note that the perfect fluid modes of stratifled mod- 
els generally have 5/3 ^ 0, because adiabatic oscillations 
generally drive fluid elements out of equilibrium, only in 
the non-stratifled case (model I) is 5/3 = strictly satis- 
fied. 

The absence/presence of g-modes in super- 
fluid/normal fluid models might seem somewhat 
surprising but one can get a better intuitive under- 
standing by considering the physical origin of these 
g-modes: a radially displaced fluid element will remain 
close to mechanical (pressure) equilibrium with its 
surroundings, but its respective values of pn and pc will 
generally differ from the surroundings (i.e. when x'^ ^ 0) 
and therefore (via the equation of state) its total density 
will also be different, resulting in a buoyant restoring 
force and a corresponding oscillation mode (in unstable 
models this restoring force will actually drive the fluid 
element still further away from its initial position, leading 



to convection). In the simple (cold) superfluid models 
considered here, each fluid element of either fluid (n or 
c) is only characterised by a single quantity, namely pn 
or pc- Displacing an element of fluid n, say, will therefore 
result not only in mechanical equilibrium (/I"), but also 
in buoyant equilibrium. This can by seen by expressing 
its density at the new position as p„ = pn{Jl^,Pc)- The 
fluid c was not displaced, therefore not only /I" but also 
Pc of the fluid element are identical to the background 
values, and so is pn- If we had allowed for an additional 
comoving quantity like entropy s, we would expect to 
find g-modes driven by a stratification in s/pn- 

It is intriguing to sec that the absence of the g-modes 
in superfluid models is accompanied by an apparent dou- 
bling of acoustic modes, but it is not obvious to establish 
a link between these different classes of modes as we are 
currently not aware of a continuous transition from a two- 
fluid to a one-fluid model (either the two fluids are locked 
together or they are not). 
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6. Conclusions 

In this paper we have tried to clarify the qualitative prop- 
erties of the eigenniode spectrum of superfluid neutron 
stars, using a simple two-fluid model. We have shown the 
important — previously somewhat overlooked — role of 
stratification for these modes. The picture has been found 
to be more complex than previous studies have suggested, 
and some of the earlier conclusions have been shown to 
apply only for non-stratified models. In particular, one 
can not generally talk about two distinct families of "su- 
perfluid" and "ordinary" modes. The system of equations 
describing two-fluid modes can not be separated in the 
case of stratified stars, and its solutions have no direct 
correspondence to the eigenmodes of the one-fluid sys- 
tem. The two-fluid modes are generally neither co- nor 
counter-moving, rather all of them are characterised by 
non-zero amplitudes of relative velocity (5 A, deviation of 
chemical equilibrium 6(3 and total density perturbation 6p. 
Also the order of the mode does not necessarily correspond 
to the number of radial nodes (as seen in Fig. ^), which is 
possible because the system is not of the Sturm-Liouville 
type. We have further confirmed earlier findings about the 
absence of g-modes in these superfluid models (Lee 1995; 
Andersson & Comer 2001a), as well as the appearance of 



avoided crossings between mode frequencies when chang- 
ing the entrainment parameter (Andersson et al. 2002). 

Given the radical difference and richer structure of the 
oscillations of superfluid neutron star models as compared 
to the simple perfect fluid models, we think that much fu- 
ture effort is needed to further clarify these properties and 
evaluate possibly observable consequences. The respective 
absence and presence of g-modes in these two different 
models is a striking example of such a potentially observ- 
able indicator of superfluidity in neutron stars. However, 
many more physical effects have to be taken into account 
in order to achieve a more realistic description of super- 
fluid neutron stars, namely the inclusion of vortex-forces 
and beta reactions, both of which will lead to dissipation. 
Furthermore, an "envelope" or an elastic crust should be 
included, and maybe most importantly, the effects of ro- 
tation and magnetic field, which add new restoring forces 
and result in a much richer spectrum of modes. Eventually, 
for a realistic study of oscillations of superfluid neutrons 
stars, one needs to work in a generally relativistic frame- 
work, as pioneered by Comer et al. (1999) and Andersson 
et al. (2002). This step is crucially important also for the 
assessment of the gravitational radiation emitted by these 
modes, and their stability /instability via the CFS mecha- 
nism (Friedman & Schutz 1978). 
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